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Abstract 

We investigate entanglement for a composite closed system endowed with a scaling property 
allowing to keep the dynamics invariant while the effective Planck constant /igff of the system is 
varied. Entanglement increases as /igff ~^ 0. Moreover for sufficiently low /igff the evolution of the 
quantum correlations, encapsulated for example in the quantum discord, can be obtained from the 
mutual information of the corresponding classical system. We show this behavior is due to the 
local suppression of path interferences in the interaction that generates the entanglement. This 
behavior should be generic for quantum systems in the classical limit. 
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Entanglement is a distinctive feature of quantum mechanics, "i/ie one that enforces its 
entire departure from classical lines of thought [1]. Its understanding has tremendously 
progressed in the last decade, due essentially to a vast amount of work regarding the con- 
struction and properties of entangled qubits in view of possible applications in quantum 
information [2^. In a more general context, any dynamical interaction between quantum 
particles leads to entanglement, that stands as a formidable obstacle to account for the 
emergence of classical word. Explaining the unobservability of entanglement in the classical 
limit is one of the aims of the decoherence programme (3| . 

Somewhat more modestly, several recent works [4] have studied in semiclassical systems 
the link between the generation of entanglement and the dynamics of the corresponding 
classical system, including in experimental realizations (5I]. The numerical and analytical 
results obtained so far indicate that the entanglement dynamics in quantum systems having 
a classically chaotic counterpart sharply differs from those whose classical counterpart is 
regular, though this difference is dependent on the specificities of the considered systems 
(types and strengths of the coupling, choice of initial states, etc.). It has been argued 
that a proper understanding of the connection between the classical dynamical regime and 
entanglement hinges on employing systems in which the same physical process generates the 
dynamics in the classical system and entanglement in its quantum counterpart. 

An intriguing question studied in this paper concerns the behavior of entanglement in 
these systems when the typical actions of the system grow with respect to h. Then the size 
of the Hilbert space increases and the quantum-classical correspondence improves. Moreover 
if the system dynamics can be kept invariant while the actions increase, an effective Planck 
constant /igfr can be defined and entanglement can be studied as f^eff ^ 0. We will see that 
entanglement indeed increases with the size of the Hilbert space in agreement with previous 
findings on entangled Bose-Einstein condensates [7] . Maybe more surprisingly for sufficiently 
low f^eff the evolution of the entanglement measure is given by probabilities obtained from the 
classical dynamical evolution, irrespective of the dynamical regime. A consequence discussed 
below is whether in the fieff ~^ limit the quantum information encoded in the pure state 
density matrix becomes indiscernible from the classical information contained in a mixed 
density matrix yielding the same reduced dynamics. 

Let us consider bipartite entanglement generated by repeated inelastic scattering of two 
particles. To set the model, let us assume a light structureless particle and a heavy rotating 
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FIG. 1: (a) The entanglement rate as a function of the number of kicks is given for quantum 
systems characterized by A: = 0.25 and /leff — 0.5, 0.25, 0.1, 0.01 (from bottom to top). The dashed 
hues result from applying the scaling relation ([4]) to the 3 upper curves, (b) Top: The black ring 
encircling the N axis is the initial classical distribution of J with J = 10 (corresponding to the 
quantum case with /leff — 0.1), centered on A/'o = T = 50. Bottom: The distribution at t = 10 
kicks (black dots), showing a spread. The surface of section for k — 0.25 is shown in yellow (light 
gray). 

particle, modeled by a symmetric top with angular momentum and energy = N{N + 
l)/2/, / denoting the moment of inertia. The scattering potential is taken to be a contact 
interaction so that the light incoming particle receives a kick when it hits the rotating top. 
The conservation of the total angular momentum T = N + J where J is the light particle 
angular momentum, imposes that after the collision the rotating top is left with an angular 
momentum A^' obeying T — J < N' <T -\- J where we have assumed J' = J. The probability 
of the transition N ^ is obtained from the scattering matrix elements iS^Nff . Finally to 
account for repeated scattering we need an attractive long-range field between both particles: 
we will assume the particles have opposite electric charge. Note that this model can be seen 
as a two-particle extension of the standard kicked top well-known in quantum chaos [sl, . 

Starting from an initial product state \iJjo) = |Fo"(6o)) |A^o) where |i^o~) depicts an in- 
coming wavepacket of the light particle with mean energy eg traveling towards the rotating 
top in state |A^o) , entanglement is generated as soon as the first collision takes place. The 
outgoing wavefunction is then given by the superposition Snnq [^^(^a^)) |^) where the 
dependence of e on is due to the conservation of energy, = E — E^^ with = eg + E^q 
being the total energy. The scattered wavepackets are later turned back by the attractive 
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field and are treated as newly incoming waves. The pure state density matrix 

p{t) = U{t,to)\iJo){HU+it,to) (1) 

is obtained by writing the evolution operator U in terms of the scattering eigenstates of the 
Hamiltonian 

= J2ZNiE)\F-{er,))\N) + J2 

N NN' 

Z-{E)Sr,,r,\F+{er,'))\N') (2) 

where are coefficients obtained by applying the asymptotic boundary conditions. The 
maximal number of entangled states is given by the number of scattering channels n = 2 
The amount of entanglement will be estimated through the entropy of the reduced density 
matrix. We will employ the linearized form 

H{t) = -^{l-Tv^p%{t)) (3) 

that becomes a good approximation for large n. pN{t) is the reduced density matrix obtained 
by tracing over the light particle's degrees of freedom. Note H = 1 for a maximally entangled 
state. For convenience we set from now on to = ^e/"^ where is the period of the mean 
energy orbit; then the collision times are t = qr^ with q being an integer. 

The classical version of the model can be formally obtained by employing the semiclassical 
link jsl between the deflection angle (f) produced by the torsional motion and the eigen- 
phaseshifts 6 of the S^-matrix: in the top's reference frame, each kick rotates J by an angle 
(j) = kJ±/J = d5/dJ± where J± is the projection of J on the unit axis Nj_ perpendicular to 
N. k is the strength of the kick; a given k corresponds, via the semiclassical relation, to a 
given S^-matrix, i.e. Snn' = SNN'{k)- The classical orbit of the light particle between two 
scattering events induces a rotation of J around N by an angle 2tit^/tn where tn is the top 
rotation period. A surface of section is obtained by plotting the position of J after each kick 
(see Figs 1(b) and 3). The crucial observation is that the surface of section only depends 
on k and on r^jr^'. J and T (which are action variables) can be increased at will, say 
by division by f^eff, but the dynamical map stays constant provided E and / are adjusted 
accordingly. For a long-range central field this also implies dividing the radial action Wr 
of the light particle by the same constant H^f^. Hence multiplying N^J^T and Wr by the 
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FIG. 2: Each panel shows the entanglement evolution as computed from the entropy H{t) (top 
curve), along with the mutual information M{t) of the classical counterpart (upside-down curve), 
where t is computed at the discrete kick times t = qr^. The parameters are J = 100, T = 500 
(corresponding to /igff — 0.01 for the quantum system) and (a): A: = 1, A/'o = 402, (b): /c = 1, A^o = 
430, (c): k = 0.01, No = 498, (d): k = 10, Nq = 460. 

common factor l//^eff 1 is tantamount to studying the limit h ^ without modifying the 
underlying dynamics. Note that the number of entangled states n also scales with l/f^eff- 

Fig. 1(a) displays the entanglement evolution for the quantum two-particle kicked top 
with k = 0.25 for different values of H^f^ (we employ atomic units and set h = 1). The light 
particle's initial distribution is a Gaussian wavepacket localized far from the symmetric top 
with its mean initial momentum directed towards it. The entanglement increases dramati- 
cally as decreases, despite the fact that the initial state takes a smaller relative area on 
the sphere. To first order, this is a consequence of the similarity transformation: on the one 
hand pN{t) is by definition a convex combination of projectors |A^) (A^|, and on the other 
hand in the semiclassical approximation the projection of pN{t) on the unit sphere (at kick 
times t = qr^) covers the same area irrespective of Hq^. Let m be the number of projectors 
|A^) (A^l (out of total of n) projecting in this area for some H^f^ and that number for 
another choice of fi^^^ < h^^. Then m/n = m! jrd from which it follows that for situations 
corresponding to the maximal entanglement (uniform distribution in that region) there is a 
simple scaling relation for the purity 1 — H yielding 



H'it) 



^eff 



(4) 



As expected entanglement increases with the number of available quantum states. 

The classical evolution analogue of the quantum problem leading to Fig. 1(a) is obtained 



by taking an initial uniform distribution of angular momentum Nq — ^ < N < Nq + 
i.e. we cut the sphere along the N axis into n slices of width 1/n. Each slice is centered so 
that the projection Jy on the Ny axis matches integer values of A^. Thus the initial classical 
distribution corresponding to |A^o) (^oL is represented on the sphere by the ring centered at 
J||(A^o) = T-N|| — A^o- The light particle initial distribution is the same Gaussian employed in 
the quantum problem (the role of this distribution is to give a statistical weight, depending 
on the initial energy of the light particle, to each J lying in the initial ring). We compute 
numerically the evolution (torsion during the kick and rotation during the orbital excursion) 
for each J of the initial distribution from kick to kick. The classical probability pfi{t = qr^) of 
finding the top with an angular momentum in the interval A^^ = ^~^^] ^^er q kicks 

is obtained by counting the vectors J whose projection falls in the corresponding interval. 
The probabilities after the first kick are given by the transition probability P^^A^^ A^). 
The classical probabilities after q kicks obey the recurrence relation 

PNiQT,) = E P^'^^^' ^ ^^)Pn' [(<! - 1) ^e] • (5) 
From these probabilities one can define the quantity 

Miqn) = ^(^l-J2 [PNiQre)]'^ (6) 

that can be understood equivalently as the linear entropy for the total system or as the 
(linearized) classical mutual information fl2|, quantifying the amount of mixing among the 
n slices induced by the kicks. 

Fig. 2(a)-(d) shows in the top panel the entanglement evolution as given by H{t = qr^) for 
a choice of system parameters (coupling k and initial state) all corresponding to f^eff ^ 1/100, 
two orders of magnitude smaller than the hard quantum case f^eff ^ 1 (which is the typical 
value for qubits), but still considerably larger than typical values characterizing classical 
actions. The bottom panel in each plot shows M[qT^) obtained from the classical probabil- 
ities through Eq. ([6]). The good agreement between H{t) and the time-dependent classical 
probabilities holds for classically chaotic and regular regimes alike as can be inferred from 
Fig. 3, displaying the corresponding surfaces of section along with the classical distributions 
whose spread along the Ny axis accounts for the entanglement evolution. 

Although quantifying entanglement by means of classical probabilities might appear sur- 
prising at first sight, we expect this behavior to be generic for semiclassical systems that 
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undergo a loss of phase coherence. This is indeed the first ingredient by which the classical 
M[qT^) can account for H{t). The second ingredient is the semiclassical approximation itself 
that allows to express operator matrix elements in terms of classical quantities (the action 
and the density of paths). For the system under consideration we start by writing Eq. ([2j) 
in the form \il){E)) = J2n ^n{E) \F{eN)) \N) where |F) is a standing wave obtained by 
combining the |F^) and Bn{E) = J^n^ SNN'Z-,{E)e'^^^^^'^^-^^/^;WP'' is the radial action 
of the classical periodic orbit in the attractive field. Then Eq. ([T]) takes the form pi 



p{t = qr,) = \N) {N'\ e'^^^^-^^'^' 



NN' 



/3^(t)/3;^,(t)|F(6^))(F(6^0l 



(7) 



with 



NN' 



N' 



-iejsft 



Z-,{E) (^(^)|^o 



^i{Wr(eN)-^)/2 



(8) 



(keep in mind cn = E — En when taking the sum). The reduced density matrix is readily 
derived as pN{t = grj = |A^) {N\pN{t) with 

2 



PN{t = qr^ 



SNN'CN'i'^) 



N' 



(9) 



C^/(t) is defined by the term between [..] in Eq. ([8]) (i.e. by excluding the phase term in 
the sum). |C7v^(t)|^ represents the probability on the incoming channel A^' just before the 
collision, whereas ASATAr^CivM weight of the outgoing wave right after the collision 

{t = qr^)] in the semiclassical limit this is the same as the weight {f^N^ at the apogee half 
a period later [14]. It follows that (^{qr^) = Pn [{q — 1) ^e] • Finally we recall \^ that in the 
semiclassical regime the S'-matrix elements are given to first order in H by 



Sj^j^, = ANN'e'^^^'/^ with \ANN'f - "^"^ 



P^'^An' A 



N) 



(10) 



and Snn' is the classical action (the boundary conditions for the conjugate momenta obey 
Po{t — oc) A^' po{t +oc) A^). As SNN'/fi oo the phase terms expi{SNNi 
Snn2)/^ Eq. ([9j) oscillate wildly, while the amplitudes Ann' are of the same order of 
magnitude. As a result these off-diagonal terms are suppressed, so we may keep only the 
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FIG. 3: Classical evolution and surfaces of section for the parameters corresponding to the plots of 
Fig. 2. (a)-(c): for A: = 1 (mixed phase-space), (a) gives the initial distributions for Nq = 402 (blue 
[grey] ring) and A^o = 430 (black ring); (b) displays the distributions after 10 kicks, (c) at t = 25. 

(d) : position of the initial distribution defined by A^o = 498 for k = 0.01 (regular dynamics). 

(e) -(f): Distributions for A; = 10 (mostly chaotic phase-space) and A^o — 460 at t = [(e)] and at 
t = 5 kicks [(f)]. 

terms with Ni = N2. Eq. ([9]) becomes 



Comparing with Eq. ([5]) and given that the initial conditions are identical in the quantum 
and classical problems, we see that provided the approximations employed hold, the en- 
tanglement entropy H{t = qr^) becomes identical to the mutual information M{qr^) of the 
corresponding classical system given by Eq. ([6]), thereby explaining the numerical results 
displayed in Fig. 2. 

A remarkable consequence of the present results concerns the classical values taken by 
quantifiers of quantum correlations. For example the quantum discord D{p) [l^ widely 
employed in the context of qubit density matrices, measures the quantum information that 
can only be extracted by joint measurements on both subsystems. D{p) vanishes if the 
state has only classical correlations. Here D{p) is simply given by H{t)^ hence by the 
classical mutual information M{t). Put diflFerently, the quantum information contained in 
the entangled state - which would be the information gained by an observer making a 
measurement (for example measuring the light particle's energy e^^ projects the top to 
the rotational state \N^)) - is given by the ignorance spread arising from the dynamical 
evolution of the corresponding classical system. 

The conjunction of ubiquitous entanglement as the classical limit is approached, and the 
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role played by classical probabilities to account for quantum correlations in that limit allows 
us to speculate whether from a statistical perspective, entanglement may be converted within 
the closed system into classical correlations, without the need to invoke couplings with an 
additional system (e.g. an environment). Indeed p{t) cannot operationally be distinguished 
from the density matrix p^^{t) = J2Pn\^) {^\\-^{^n)) {F{eN)\ containing only classical 
correlations (and for which D{p^^) = 0): the reduced density matrices obtained from p and 
p^^ are identical, and as ^ the coherences (in the "pointer basis" \ F{eN)) \N)) of typical 
two-particle observables would lead to interference patterns with vanishing (and therefore 
undetectable) wavelengths 

To sum up, we have investigated entanglement evolution when the entanglement is gener- 
ated by a dynamical localized interaction in a quantum system having a well defined classical 
counterpart. We have seen that entanglement increases, irrespective of whether the under- 
lying dynamics is regular or chaotic, as typical actions grow relative to h. The quantum 
correlations are then given by the mutual information of the corresponding classical system. 
The present results could contribute to a better understanding of the role played by quantum 
information in the classical limit. 
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